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Abstract 


Using a renormalization method, we study the critical behavior for intermit¬ 
tency in two coupled one-dimensional (ID) maps. We find two fixed maps of 
the renormalization transformation. They all have common relevant eigen¬ 
values associated with scaling of the control parameter of the uncoupled ID 
map. However, the relevant “coupling eigenvalues” associated with coupling 
perturbations vary depending on the fixed maps. It is also found that the two 
fixed maps are associated with the critical behavior in the vicinity of a critical 
line segment. One fixed map with no relevant coupling eigenvalues governs 
the critical behavior at interior points of the critical line segment, while the 
other one with relevant coupling eigenvalues governs the critical behavior at 
both ends. The results of the two coupled ID maps are also extended to many 
globally-coupled ID maps, in which each ID map is coupled to all the other 
ones with equal strength. 
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I. INTRODUCTION 


An intermittent transition to chaos in the ID map occurs in the vicinity of a saddle- 
node bifurcation [Q], Intermittency just preceding a saddle-node bifurcation to a periodic 
attractor is characterized by the occurrence of intermittent alternations between laminar 
and turbulent behaviors. Scaling relations for the average duration of laminar behavior 
in the presence of noise have been first established |@] by considering a Langevin equation 
describing the map near the intermittency threshold and using Fokker-Plank techniques. 
The same scaling results have been later found @ by using the same renormalization-group 
equation for period doubling with a mere change of boundary conditions appropriate to 
a saddle-node bifurcation. 

Recently, efforts have been made to generalize the scaling results of period doubling 
for the ID map to coupled ID maps P-[TI[|, which are used to simulate spatially extended 


systems with effectively many degrees of freedom 


It has been found that the critical 


scaling behaviors of period doubling for the coupled ID maps are much richer than those 
for the uncoupled ID map HI- These results for the abstract system of the coupled 
ID maps are also confirmed in the real system of the coupled oscillators [^. In a similar 
way, the scaling results of the higher period p-tuplings (p =,3,4,...) in the ID map are 
also generalized to the coupled ID maps |^. Here we are interested in another route to 
chaos via intermittency in coupled ID maps. Using a renormalization method, we extend 
the scaling results of intermittency for the ID map to coupled ID maps. 

This paper is organized as follows. In Sec. || we introduce two coupled ID maps and 
discuss their symmetry. Bifurcations associated with stability of periodic orbits are also 
discussed there. In Sec. we employ the same renormalization method |]^ developed for 


period doubling in coupled ID maps and study the critical behavior for intermittency in 
two coupled ID maps. We thus find two fixed maps of the renormalization transformation. 
They have the relevant eigenvalues associated with scaling of the control parameter of the 
uncoupled ID map as common ones. However, the relevant “coupling eigenvalues” (CE’s) 
associated with coupling perturbations vary depending on the fixed maps. These two fixed 
maps are also found to be associated with the critical behavior near a critical line segment in 
Sec. |I^. One fixed map with no relevant CE’s governs the critical behavior at interior points 
of the critical line, while the other one with relevant CE’s governs the critical behavior at 
both end points. We also study the critical behavior for intermittency in many globally- 
coupled ID maps in Sec. 0. Globally-coupled systems, in which each element is coupled to 
all the other ones with equal strength, appear naturally in broad branches of science ||I5[ . 
The results of two coupled maps are extended to this kind of many globally-coupled maps. 
Finally, a summary is given in Sec. B 


II. TWO COUPLED ID MAPS 

After briefly reviewing the intermittency in case of the ID map, we introduce two coupled 
ID maps and discuss their symmetry. Bifurcations associated with stability of periodic orbits 
are also discussed. 

We first recapitulate the intermittent transition to chaos JIJ-Q in a ID map with one 
control parameter A, Xt+i = u{Xt) (t denotes a discrete time). A pair of orbits with period p 
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appears via saddle-node bifurcation, as the control parameter A exceeds a threshold value Ac. 
One periodic orbit is a stable attractor, while the other one is an unstable repeller. However, 
as A decreases below Ac, the two periodic orbits disappear, and an intermittent chaotic 
attractor, characterized by the occurrence of intermittent alternations between laminar and 
turbulent behaviors, appears. 

One can easily explain the intermittency geometrically as follows. The curve of the 
equation Y = u^^\X) is the pth iterate of u] has new 2p intersection points with the 
Y = X line for A > Ac, which collapse into p points tangent to the Y = X line for A = Ac 
(i.e., we have p hxed points, X^ = u^^\Xl), t = 1, ... ,p ioi A = Ac). However, as A 
decreases below Ac the curve no longer touches the Y = X line so that a “channel” appears 
in the immediate vicinity of each point tangent to the Y = X line at A = Ac. If one orbit 
point falls close to the entrance to one of the channels under repeated iterations of it 
would take many iterations to go through the channel, which corresponds to the “laminar 
phase.” After slow passage through the channel, the iterates of move wildly until they 
return to one of the channels. This corresponds to the “turbulent phase.” Thus the laminar 
and turbulent phases appear intermittently. 

The nearer A is to Ac, the longer the averaged laminar time. To obtain the average 
duration of the laminar phase, consider the pth iterate u^^\X) in the immediate vicinity of 
one of the channels. Shifting the origin of coordinate X to a hxed point X* of u^^\X) for 
A = Ac, we have 

xt+i = f{xt) = u^^\xt + X*) - X* ^ xt + alxtl"" + e, z > 1, ( 2 . 1 ) 


where x = X — X*, a is a constant, and e is a control parameter proportional to A — Ac. 
Using a map / of the form (p.l|) , it has been found in that the mean duration of the 
laminar phase l{e) varies as Hence the tangency-order z determines the universality 

classes, because I depends on the order z. In this paper, we consider only the analytic case 
of even order z (z = 2,4, 6,...). 

We now consider a map M consisting of two identical ID maps coupled symmetrically. 


M : 


W+i = W{Xt, Yt) = u{Xt) + v{Xt, Yt), 
Ui+i = WiYt,Xt) = u{Yt)+v{Yt,Xt), 


( 2 , 2 ) 


where n is a coupling function obeying a condition. 


n(X,X) = 0 for any X (2.3) 

The map ( |2.2| ) is called a symmetric map because it has an exchange symmetry such 
that 


cT-Wa(Z) = M(Z) for all Z, 


(2.4) 


where Z = {X,Y), a is an exchange operator acting on Z such that o'(Z) = {Y,X), and 
a~^ is its inverse. The set of all hxed points of a forms a synchronization line U = X in 
the state space. It follows from Eq. (|2.4| ) that the exchange operator a commutes with the 
symmetric map M, i.e., aM = Ma. Thus the synchronization line becomes invariant under 
M, i.e., if a point Z lies on the synchronization line, then its image M{Z) also lies on it. An 
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orbit is called a(n) (in-phase) synchronous orbit if it lies on the synchronization line, i.e., it 
satisfies 


Xt = Yt = X; for all t. 


(2.5) 


Otherwise, it is called an (out-of-phase) asynchronous orbit. Here we study the intermittency 
associated with a saddle-node bifurcation to a pair of synchronous periodic orbits, which 
can be easily found from the uncoupled ID map, X^j^^ = m(Xj*), because of the coupling 
condition (|2.3|) . 

Stability of a synchronous orbit of period p is determined from the Jacobian matrix J of 
(pth iterate of M), which is given by the p product of the linearized map DM of the 
map (p.2|) along the orbit 


t=l 

^ -n ( - v(.x;) v(x;) \ 

M1 v(x;) u'{x;) - v(x;) j • 


(2.6) 


where u'{x) = du{X)/dX and V{X) = dv{X,Y)/dY\Y=x', hereafter V{X) will be referred 
to as the “reduced coupling function” of g{x,y). The eigenvalues of J, called the stability 
multipliers of the orbit, are gievn by 

x, = flu\x:), \^ = f{\u'(Xl)-2V{X:)]. (2.7) 

t t 

Note that Ai is just the stability multiplier for the case of the uncoupled ID map and the 
coupling affects only A 2 . 

A synchronous periodic orbit is stable when both multipliers he inside the unit circle, 
i.e., |Aj| < 1 for j = 1 and 2. Thus its stable region in the parameter plane is bounded by 
four bifurcation lines, i.e., those curves determined by the equations \j = ±1 (j = 1,2). 
When a multiplier Aj increases through 1, the stable synchronous periodic orbit loses its 
stability via saddle-node or pitchfork bifurcation. On the other hand, when a multiplier 
Xj decreases through —1, it becomes unstable via period-doubling bifurcation. (For more 
details on bifurcations, refer to Ref. [pR| .) 


III. RENORMALIZATION ANALYSIS OF TWO COUPLED MAPS 


Here we are interested in intermittency just preceding a saddle-node bifurcation. The 
intermittent transition to chaos occurs near a critical line segment, as will be seen in Sec. 
Using the same renormalization method [pi]| , p!T| for period doubling with a mere change of 
boundary conditions appropriate for a saddle-node bifurcation, we generalize the ID scaling 
results for intermittency to the case of two coupled ID maps. We thus find two hxed maps of 
the renormalization transformation and obtain their relevant eigenvalues. One hxed map is 
associated with the critical behavior inside the critical line, while the other one is associated 
with the critical behavior at both ends. 
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To study the intermittent transition to chaos near a saddle-node bifurcation to a pair of 
synchronous orbits of period p, consider the pth iterate of M of Eq. (p.2|) , 

= W^P\Yt,Xt), (3.1) 

where satishes a recurrence relation 

W^P\X,Y) = W{W^p-^\X,Y),W^p-^\Y,X)). (3.2) 

The function can be decomposed into the uncoupled part u^p'* and the remaining 

coupling part, i.e., 

W^P\X, Y) = u^P\X) + [W^P\X, Y) - m(^’)(X)]. (3.3) 

When the control parameter A of the uncoupled ID map is equal to the threshold value 
Ac for the synchronous saddle-node bifurcation, we have p synchronous hxed points of M^p^ 
such that Y^* = X^ = u^p\X^) for f = 1,... ,p. Shifting the origin of coordinates (X, Y) to 
one of the p synchronous hxed points (X*, Y*) {Y* = X* = u^p\X*) for A = Ac), we have 

rp . f xt+i = F{xt,yt) = f{xt)+g{xt,yt), 

■ \ yt+i = F{yt,Xt) = fiyt) + giyt,Xt), 

where 

f{x) = u^P\x + X*) - X*, 
g{x,y) = W^P\x + X%y + Y*)-u^\x + X*). 

Since a synchronous saddle-node bifurcation occurs at the origin (0, 0) for A = Ac in case of 
the map T, the uncoupled part / for the critical case A = Ac satishes 


(3.4) 

(3.5) 

(3.6) 


/(0) = 0 and /'(0) = 1. (3.7) 

Note also that the coupling function g{x,y) satishes the coupling condition (|2.3|) , i.e., 

g{x,x) = 0 for any X. (3.8) 


We employ the same renormalization transformation |]T0|,pT| as in the period-doubling 


case with changed boundary conditions (|3.7|) . The renormalization transformation Af for a 
coupled map T consists of squaring (T^) and rescaling (B) operators: 


A/'(T) = BT^B-^ 


(3.9) 


Since we consider only synchronous orbits, the rescaling operator is of the form. 


B = 




(3.10) 


where a is a rescaling factor. 

Applying the renormalization operator Af to the coupled map (|3.4|) n times, we obtain 
the u-times renormalized map of the form. 
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T 

-L n. 


Xt+l = Fn{xt,yt) = fn{xt) +gn{xt,yt), 
yt+l = FniyuXt) = fniyt) + gn{yt,Xt). 


(3.11) 


Here /„ and gn are the uncoupled and coupling parts of the n-times renormalized function 

(3.12) 


Fn, respectively. They satisfy the following recurrence equations |T0|JTl 


,x, 


fn+l{x) = afnifni-)), 

a 

gn+l{x,y) = afnifni-) + gni-, “)) 
a a a 

+agn{fn{-)+9n{-,-)Jn{-) 
a a a a 

,y a;,. „ , „ ,a;.. 

+9n{-,-)) -afnifni-))- 
OL a a 


(3.13) 


Then Eqs. ( p.l2|) and ( |3.13| ) define a renormalization operator F of transforming a pair of 
functions, {f,9)', 


fn+l 

9n+l 


= n 


fn 

9n 


(3.14) 


A critical map with the nonlinearity and coupling parameters set to their critical 
values is attracted to a fixed map T* under iterations of the renormalization transformation 

AT, 




xt+i = F*{xt,yt) = f*{xt) + 9* {xt,yt), 
yt+i = F*{yt,xt) = r{yt)+9*{yuxt). 


(3.15) 


Here {f*,9*) is a hxed point of the renormalization operator TZ, which satisfies {f*,9*) = 
9*)- Note that the equation for f* is just the hxed-point equation for the intermittency 
with boundary conditions (|3.7|) in the uncoupled ID map. It has been found in p] that 


r{x) = (r[l -{z- l)a(r"-']-^A-i) 


= X + ax + -a X 


2„2z-l 


+ ■ ■ ■ (a : arbitrary constant) 
is a fixed point of the transformation (|3.12|) with 


(3.16) 


(3.17) 


[As mentioned in Sec. ||, we consider only the analytic case of even order z {z = ...).] 

Consequently, only the equation for the coupling hxed function g* is left to be solved. 

However, it is not easy to directly solve the equation for the coupling hxed function. We 
therefore introduce a tractable recurrence equation for a reduced coupling function of the 
coupling function g{x,y), dehned by 


G(.) ^ 


dy 


(3.18) 


y=x 


Diherentiating the recurrence equation (|3.13|) for g with respect to y and setting y = x, we 
have 
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(3.19) 


Gn+i(a:) 


[/'(/n(-))-2G„(A(-))]G.(-) 

a a a 

+Gnifni-))fLi-)- 

a OL 


Then Eqs. ( |3.12|) and ( p.l9| ) define a “reduced renormalization operator” 7?. of transforming 
a pair of functions (/, G) |1TmpT| : 


f fn+l 

\ Gn+l 



(3.20) 


We now look for a fixed point of Iz, which satisfies = iZ{f*,G*). Here 

G* is just the reduced coupling fixed function of g* (i.e., G*{x) = dg*{x,y)ldy\y=x)- Using 
a series-expansion method, we find two solutions for G*\ 

G;(i) = i|l + zax‘-' + 0(2 - + ...], (3.21) 

G],(x) = + ...], (3.22) 


where a and b are arbitrary constants. Here we are able to sum the series in Eq. (|3.21|) and 
obtain a closed-form solution, 


GKx) 



(3.23) 


However, unfortunately we cannot sum the series in Eq. ( |3.22|) , except for the cases 6 = 0 
and z where we obtain closed-form solutions. 


Ghix) 


0 for 6 = 0, 
^[f*'ix) — 1] for b = z. 


(3.24) 


Consider an infinitesimal perturbation {h, $) to a fixed point (/*, G*) of the reduced 
renormalization operator TZ. We then examine the evolution of a pair of functions (/* -|- 
h, G*-l-<h) under TZ. Linearizing TZ at its fixed point, we obtain a reduced linearized operator 
C transforming a pair of perturbations (h, $): 



where 

hn+i{x) = [Cihn]{x) = (3.26) 

a a a 

tn+l)!) = li:2tn](2^) + [G'iK][x), (3.27) 

[£2<1>„](i)= [/•'(/•(-)) - 2G*(/*(-))]<l>„(-) + [/•'(-) - 2G-(b]1>„(/*(-)), (3.28) 

a a a a a a 

[c,K]{x)=[r'\f\-))G\-)-2G*\f\-))G\-) + G*\r{-))r\^K{-) 

a a a a a a a 

+G%-mn-))+G*in-)m-). ( 3 . 29 ) 

a a a a 


7 











Here the prime denotes a derivative. Note that, although hn couples to both hn+i and <hn+i, 
couples only to ‘h^+i- From this reducibility of C into a semiblock form, it follows that 
to obtain the eigenvalues of C it is sufficient to solve the eigenvalue problems for Ci and £2 
separately. The eigenvalues of both £1 and £2 give the whole spectrum of £. 

A pair of perturbations {h*, <h*) is called an eigenperturbation with eigenvalue A, if 



that is, 

Xh*{x) = [£ih*](a;). 


(3.30) 


(3.31) 

(3.32) 


We hrst solve Eq. ( p.31|) to hud eigenvalues of £ 1 . Note that this is just the eigenvalue 
equation for the ID map case. The complete spectrum of eigenvalues and the corresponding 
eigenfunctions have been found in Refs. [^. The form of the eigenvalues is A„ = 

(n = 0,1, 2,...). Hence the hrst z eigenvalues with n < z are relevant ones. The marginal 
one A,j is associated with the arbitrary constant a in f*{x), and all the other ones with n > z 
are irrelevant. Although the eigenvalues A„’s of £1 are also eigenvalues of £ as mentioned 
in the preceding paragraph, [h*, 0 ) itself cannot be an eigenperturbation of £ unless £3 is a 
null operator. 

We next consider a perturbation of the form (0, <h) having only the coupling part. If a 
reduced coupling perturbation satishes 


A<l>*(a;) = [£ 2 <h*](x) 

= [/*'(/*(-))-2G*(r(-))]$*(-) 

a a a 

+ [r'(-)-2G*(-)]r(r(-)), (3.33) 

a a a 


then it is called a reduced coupling eigenperturbation with coupling eigenvalue (CE) A. In 
this case (0, d**) becomes an eigenperturbation of £ with CE A. 

We hrst consider the case of the solution G*j{x) = A/*'(a:). In this case the linearized 
operator £2 of Eq. (|3.33|) becomes a null operator because the right-hand side of the equation 
becomes zero. Hence there exist no relevant CE’s associated with coupling perturbations, 
and consequently the hxed point (/*, G}) has only relevant eigenvalues of £1 associated with 
scaling of the control parameter of the uncoupled ID map. As will be seen in Sec. |I^, the 
intermittent transition to chaos occurs near a critical line segment. The critical behaviors 
near interior points of the critical line segment are essentially the same as those for the ID 
case. In fact, a pair of critical functions (/c, Gc) at any interior point is attracted to the 
same hxed point (/*,G/) under iterations of 'JZ. Hence the hxed point (/*,G/) governs the 
critical behaviors inside the critical line segment. (For details refer to the next section.) 

Second, consider the case of the solution G}j of Eq. (|3.22|) . Using a series-expansion 
method, we hnd the complete spectrum of CE’s and the corresponding eigenfunctions. An 
eigenfunction <I)*(a:) can be expanded as follows: 


00 


1=0 


(3.34) 
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Substituting the power series of f*{x), f*'{x), G*jj{x) and <h*(x) into the eigenvalue equation 
( p.33| ), it has the structure 


\ck = ^MkiCu fc,/= 0,1,2,... . 

i 


(3.35) 


Note that each q (/ = 0,1, 2,...) in the right-hand side is involved only in the determination 
of coefficients of monomials x^ with k = I + m{z — 1) (m = 0,1,2,...). Hence A4 becomes 
a lower triangular matrix. Its eigenvalues are therefore just diagonal elements: 

= = ^ = fc = 0,l,2,.... (3.36) 

a'* 

The hrst {z — 1) eigenvalues A^’s for 0 < k < z — 2 are relevant ones. The marginal eigenvalue 
A^-i is associated with the arbitrary constant b in G*jj{x), and all the other eigenvalues for 
k > z — 1 are irrelevant. 

Each eigenfunction <hfc(a;) with CE A*, (/c = 0,1, 2,...) is of the form, 

<hfc(a;) = x’^ll + {z - b+ ^)ax"'~^ H-]. (3.37) 

In case of the largest CE Aq = 2, we are able to sum the series in Eq. (|3.37|) for 6 = 0 and ^ 
and hnd *ho(^) closed form: 




f*'{x) for 6 = 0 , 
1 for b = z. 


(3.38) 


We can also sum the series in Eq. (|3.37|) for all the irrelevant cases (i.e., the cases k > z) 
and hnd the closed-form eigenfunctions. 




[X 


= - 2 Gj^(a;)][/*'(a;) -x'], I = k - z + 1 = 1,2,... 


(3.39) 


which are associated with coordinate changes 

As we shall see in the next section m the second reduced coupling function G*jj{x) 
governs the critical behavior near both ends of a critical line segment, because a pair of 
critical functions (/c, Gc) at each end point is attracted to the hxed point {f*,G}j) under 
iterations of 77. (For details refer to the next section.) 


IV. SCALING BEHAVIOR NEAR A CRITICAL LINE 


We choose the uncoupled ID map in two coupled ID maps M of Eq. (|2.2|) as 

u{X) = 1- AX^, (4.1) 

and consider a dissipative coupling case in which the coupling function is given by 

v(X,Y) = IHY)-u(X)1 (4,2) 
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Here c is a coupling parameter. There exists a critical line segment associated with intermit- 
tency on the synchronous saddle-node bifurcation line A = Acin the c — A plane. Inside the 
critical line segment, the critical behaviors are governed by the fixed point (/*, G*j) with no 
relevant CE’s, and hence they are essentially the same as those for the ID case. On the other 
hand, the critical scalings near both ends are governed by the hxed point (/*, G*jj) with one 
relevant CE. As c passes through both ends, the coupling leads to desynchronization of the 
interacting systems. Hence the synchronous attractor ceases to be an attractor and a new 
asynchronous attractor appears outside the critical line segment. Thus a transition from a 
synchronous to an asynchronous state occurs at both ends of the critical line segment. 

As an example, we consider the saddle-node bifurcation to a pair of synchronous orbits 
with period p = 3 occurring for A = Ac = 1.75. To study the intermittency associated 
with this bifurcation, we first consider the third iterate of M [see eq. (0)], and then 
shift the origin of coordinates {X, Y) to one of the three synchronous fixed points {X*, Y*) 
for A = Ac \Y* = X* = Thus we obtain a map T of the form ( |3.4| ), where the 

uncoupled and coupling parts / and g are given by 

f{x)=u^^\x + X*)-X\ (4.3) 

g{x, y) = W^^\x + X\y + Y*)- u^^\x + X*). (4.4) 


Near the region of the synchronous saddle-node bifurcation, f{x) can be expanded about 
X = 0 and A = Ac, 


f{x) ~ X -t- ax^ -f e. 


(4.5) 


where a = f /dx‘^\x=Q,A=Ac ci^nd e = df /dA\x=Q,A=Ac (^ ~ ^c)- Hence this corresponds to 
the most common case with the tangency-order z = 2 [see Eq. (|2.1|)]. The reduced coupling 


function G{x) of g{x,y) [defined in Eq. (|3.18|) ] is also given by 


G{x) = -f{x), e = c^ 


3m -|- 3c. 


(4.6) 


Consider a pair of initial functions (/c, G) on the synchronous saddle-node bifurcation 
line A = Ac, where fc{x) is just the ID critical map and G(x) = |/c(a^)- By successive 
actions of the reduced renormalization transformation TZ of Eq. (|3.2CI|) on (/c, G), we obtain 


Ux) = a/._i(/._i(-)), G.(x) = ^ /;(x), (4.7) 

a 2 

e„= 2e„_i - e^_i, (n = 0,1, 2,...), (4.8) 


where the rescaling factor for = 2 is a = 2 , /o(x) = fc{x), Go(x) = G(x), and Cq = e. Here 
fn converge to the ID fixed function f*(x) of Eq. (|3.16|) with z = 2 as n —>• cx). 

Figure Q shows a plot of the curve determined by Eq. ([4.8|) . Two intersection points 
between this curve and the straight line e„ = e„_i are just the fixed points e* of the recurrence 
relation ( |4.8D for e. 


e* = 0, 1. (4.9) 

Stability of each hxed point e* is determined by its stability multiplier A [= (ie„/de„_i|e*]. 
The hxed point at e* = 1 is superstable (A = 0), while the other one at e* = 0 is unstable 
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(A = 2). The basin of attraction to the superstable hxed point e* = 1 is the open interval 
(0, 2). That is, any initial e inside the interval 0 < e < 2 converges to e* = 1 under successive 
iterations of the transformation ( [4.81) . The left end of the interval is just the unstable hxed 
point e* = 0, which is also the image of the right end point under the recurrence equation 
([4.8|). All the other points outside the interval diverge to the minus inhnity under iterations 


of the transformation (|4.8|) . 

It follows from the relation e = e(c) in Eq. ( |4.6|) that there exists a critical line segment 
joining two end points Q = 0 and = 2 on the synchronous saddle-node bifurcation line 
A = Ac in the c — A plane. Inside the critical line segment, any initial G{x) is attracted 
to the hrst reduced coupling hxed function G}{x) = ^f*'{x) under iterations of TZ, while 
(7(0:)’s at both ends are attracted to the second reduced coupling hxed function G*jj{x) = 0 
with 6 = 0. 

Figure ^ shows a phase diagram for the case of a dissipative coupling (|4.2|) . The diagram 
is obtained from calculation of Lyapunov exponents. For the case of a synchronous orbit, 
its two Lyapunov exponents are given by 


1 


m—1 


t=0 
T m—1 

(Jx{A,c) = ^h^ —^ ln|(l -c)M'(Xt)| = a||(A) +ln|l -c|. 
^ t=o 


(4.10) 

(4.11) 


Here cr||((T^) denotes the mean exponential rate of divergence of nearby orbits along (across) 
the synchronization line Y = X. Hereafter, a\\ (cx^) will be referred to as tangential and 
transversal Lyapunov exponents, respectively. Note also that ct|| is just the Lyapunoiv 
exponent for the ID case, and the coupling ahects only a±. 

The data points on the crj^ = 0 curve are denoted by solid circles in Fig. A synchronous 
orbit on the synchronization line becomes a synchronous attractor with < 0 inside the 

a± = 0 curve. The type of this synchronous attractor is determined according to the sign 
of (j||. A synchronous period-3 orbit with a\\ < 0 becomes a synchronous periodic attractor 
above the critical line segment, while a synchronous chaotic attractor with a\\ > 0 exists 
below the critical line segment. The periodic and chaotic parts in the phase diagram are 
denoted by P and G, respectively. There exists also a synchronous period-3 attractor with 
(T|| = 0 on the critical line segment between the two parts. 

Consider a transition to chaos near an interior point with q < c < of the critical line 
segment. Here we £x the value of the coupling parameter c and vary the control parameter 
(= Ac —A). For < 0 there exists a synchronous period-3 attractor on the synchronization 
line Y = X. However, as eA is increased from zero, the periodic attractor disappears and 
a new chaotic attractor appears on the synchronization line. As an example, see the hgure 
^(a) which shows a synchronous chaotic attractor for the case ca = 10““^. The motion on this 
synchronous chaotic attractor is characterized by the occurrence of intermittent alternations 
between laminar and turbulent behaviors on the synchronization line, as shown in Fig. |^(b). 
This is just the intermittency occurring in the uncoupled ID map, because the motion on 
the synchronization line is the same as that for the uncoupled ID case. Thus, a “ID-like” 
intermittent transition to chaos occurs near interior points of the critical line segment. 

The scaling relations of the mean duration of the laminar phase I and the tangential 
Lyapunov exponent a\\ for a synchronous chaotic attractor are obtained from the leading 
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relevant eigenvalues (= 4) of the fixed point of TZ with no relevant CE’s, like 

the ID case [§]. A map with non-zero e near an interior point of the critical line segment 
is transformed into a new map of the same form, but with a new parameter e' under a 
renormalization transformation. Here the control parameter scales as 

6^4 = 5ieA = 22e^. (4.12) 

Then the mean duration I and the tangential Lyapunov exponent a\\ satisfy the homogeneity 
relations, 

K^'a) = ^K^a), = 2a||(eA), (4.13) 

which lead to the scaling relations. 


K^a) 


-A ’ 




i^A) 


-A’ 


with exponent 


p = log 2/ log^i = 0.5. 


(4.14) 


(4.16) 


The ID-like intermittent transition to chaos ends at both ends q and Cr of the critical 
line segment. We £x the value of the control parameter A = Ac{= 1.75) and study the 
critical behaviors near the two end points by varying the coupling parameter c. Inside 
the critical line segment (q < c < c^), a synchronous period-3 attractor exists on the 
synchronization line. However, as the coupling parameter c passes through q or c^, the 
transversal Lyapunov exponent a± of the synchronous periodic orbit increases from zero, 
as shown in Fig. and hence the coupling leads to desynchronization of the interacting 
systems. Thus the synchronous period-3 orbit ceases to be an attractor outside the critical 
line segment, and a new (out-of-phase) asynchronous attractor appears. This is illustrated 
in Fig. There exist an asynchronous attractor with period 3 denoted by uptriangles 
for c = —0.0001, while an asynchronous period-6 attractor denoted by “stars” exists for 
c = 2.0001. Since the asynchronous period-3 attractor is not a symmetric one under the 
exchange operator a of Fq. (2.4), there exists also its conjugate asynchronous period-3 
attractor denoted by downtriangles. But the asynchronous period-6 attractor is a symmetric 
one under the exchange operator a. 

The critical behaviors near both ends q and Cr are the same. Since the transversal 
Lyapunov exponent a± is equal to ln|l — c| for A = Ac [see Eq. (|4.11|) ], it varies linearly 


with respect to c near both ends, (i.e., a± 


ec = Q 


c or c 


The scaling 


behavior of a_i_ is obtained from the relevant CE 62 (= 2) of the hxed point {f*,GJj) of TZ 
as follows. Consider a map with non-zero Cc, (but with ca = 0) near both ends. It is then 
transformed into a new one of the same form, but with a renormalized parameter e), under 
a renormalization transformation. Here the parameter Cc obeys a scaling law, 

e'c = (^2 Q = 2ec. (4.16) 

Then the transversal Lyapunov exponent cxj, satishes the homogeneity relation, 

ade'c) = 2ac{ec). (4.17) 
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This leads to the scaling relation, 


ffc(£c) ~ 


(4.18) 


with exponent 


V = log 2/ log52 = 1- 


(4.19) 


Like the case of interior points of the critical line segment, the scaling behavior of (T\\{eA) 
for c = Q or Cr is obtained from the relevant eigenvalue hi (= 4) of the fixed point (/*, G}j), 
and hence it also satishes the scaling relation ( 4.10 ). The scaling behaviors of the Lyapunov 
exponents ct|| and a± near both ends are thus determined from two relevant eigenvalues hi 
and h 2 of the hxed point (/*, GJj). 


V. EXTENSION TO MANY GLOBALLY-COUPLED MAPS 


Recently, much attention has been paid to dynamical systems with many nonlinear 
elements and a global coupling, in which each element is coupled to all the other elements 
with equal strength. Globally-coupled systems as the extreme limit of long-range couplings 
are seen in broad branches of science |]J]. For example, coupled nonlinear oscillators with a 
global coupling frequently occur in charge density waves [|1^, Josepson junction arrays 
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and p-n junction arrays This kind of dynamical systems with a global coupling can be 
also regarded as mean-field versions of dynamical systems with local short-range couplings. 

Here we study the critical behavior for intermittency in many-coupled ID maps with 
a global coupling. The results of two coupled maps are straightforwardly extended to the 
many globally-coupled maps. 


A. Many globally-coupled ID maps 

Consider an iV-coupled map with a periodic boundary condition: 

M : Xt+iim) = W{a^-^Xt) 

= + 1),..., X,{m - 1)), 

m = 1,..., iV, (5.1) 

where the number of constituent elements iV is a positive integer larger than or equal to 
2, Xt{m) is the state of the mth element at a lattice point m and at a discrete time t, 
X = (X(l), A(2),... ,X{N)), a is the cyclic permutation of the elements of X, i.e., aX = 
(X(2),..., X(iV), A(l)), means (m — 1) applications of a. The periodic condition 

imposes Xt{m) = Xt{rn + N) for all m. Like the two-coupled case (iV = 2), the function W 
consists of two parts: 


U-(X) = i.(A'(l)) + k(X), 


(5.2) 


where u is an uncoupled ID map and n is a coupling function obeying the condition. 
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t>(X,..., X) = 0 for any X. 


(5.3) 


Thus the iV-coupled map (El becomes: 

M : Xt+i{m) = u{Xt{m)) 

+v{Xt{m),Xt{rn + 1),..., Xtijn - 1)), 

m = 1,..., iV. (5.4) 

Here we study many-coupled ID maps with a global coupling. In the extreme long-range 
case of global coupling, the coupling function v is of the form, 

c ^ 

n(X(l),..., X(iV)) = ^ E HX{m)) - r(X(l))] for N>2, (5.5) 

where r{X) is a function of one variable. Note that each ID map is coupled to all the other 
ones with equal coupling strength c/N inversely proportional to the number of degrees of 
freedom N. Hereafter, c will be called the coupling parameter. 

The A^-coupled map M is called a symmetric map, because it has a cyclic permutation 
symmetry such that 


a-^Ma(X) = M(X) for all X, (5.6) 

where a~^ is the inverse of a. The set of all hxed points of a forms a synchronization line 
in the A^-dimensional state space, on which 

A'(l) = ... = X(N). (5.7) 

It follows from Eq. (El) that the cyclic permutation a commutes with the symmetric map 
M, i.e., aM = Ma. Hence the synchronization line becomes invariant under M. An orbit is 
called a(n) (in-phase) synchronous orbit if it lies on the synchronization line, i.e., it satishes 

Xi(l) = ■ ■ ■ = Xt{N) = X; for all t. (5.8) 

Otherwise, it is called an (out-of-phase) asynchronous orbit. Here we study the intermittency 
associated with a synchronous saddle-node bifurcation. Note also that synchronous orbits 
can be easily found from the uncoupled ID map, X^j^^ = u{X^), because of the coupling 
condition (|5.3|) . 


B. Lyapunov exponents of synchronous orbits and the critical behavior 


Consider an orbit {X^} = {Xt{m), m = 1,..., N} in many coupled maps (El)- Stability 
analysis of the orbit can be easily carried out by Fourier transforming with respect to the 
discrete space {m}. The discrete spatial Fourier transform of the orbit is: 


N 


= ^Y. e-"™’*"" A',(m) = f,(j), 

m=l 

j=0,l,...,iV-l. 


(5.9) 


14 



The Fourier transform ^t{j) satisfies (j) = — j) (* denotes complex conjugate), and 

the wavelength of a mode with index j is y for j < y and for j > y. Here ^^(0) 
corresponds to the synchronous (Fourier) mode of the orbit, while all the other with 

nonzero indices j correspond to the asynchronous (Fourier) modes. 

To determine the stability of a synchronous orbit [Xt(l) = ■ ■ ■ = Xt{N) = for 
all t], we consider an inhnitesimal perturbation {AXt{m)} to the synchronous orbit, i.e., 
Xt{m) = X^ + AXt{m) for m = Linearizing the iV-coupled map (|5.4|) at the 

synchronous orbit, we obtain: 

N 

AXt+i{m) = u'{X;)AXtim) + Y.V^^\x:)AXt{l + m - 1), (5.10) 

where 

mx) ^ v%x ). 


Hereafter, the functions will be called “reduced coupling functions” of u(X). 

Let be the Fourier transform of AXt{m), i.e., 

1 ^ 

S(,U) = rlAX,(m)] = - y 

m=l 

j = 0,l,...,iV-l. (5.12) 

Here ^^t(0) is the synchronous-mode perturbation along the synchronization line, while all 
the other ones 5^t(j)’s with no n zero indicies j are the asynchronous-mode perturbations 
across the synchronization line. Then the Fourier transform of Eq. ( |5.1U|) becomes: 

N 

l=l 

j = 0,l,...,iV-l. (5.13) 


1=1 


^^(^(m-l)x) 


1 

dv{X) 


X(1)=---=X{N)=X 


dXi 


(5.11) 


X(1)=---=X(N)=X 


Note that all the modes ^t(j)’s become decoupled for the synchronous orbit. 

The Lyapunov exponent Xj of a synchronous orbit, characterizing the average exponential 
rate of divergence of the jth-mode perturbation, is given by 

1 m—1 N 

a, = Jto -y In|,.'(A7) + yi/«(A7)e""l'-'M-'"|. (5.14) 

t=0 l=l 


If the Lyapunov exponent cxj is negative or zero, then the synchronous orbit is stable against 
the jth-mode perturbation; otherwise it is unstable. 

In case of a synchronous periodic orbit with period p, the Lyapunov exponent aj in 
Eq. (^.14| ) becomes: 
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(5.15) 


aj = - In I A,-1, 

p 

where Xj, called the stability multiplier of the synchronous orbit, is given by 


p-i 


N 




(5.16) 


t=0 


1=1 


The synchronous periodic orbit is stable against the jth-mode perturbation when Xj lies 
inside the unit circle, i.e., |Aj| < 1. When the stability multiplier Xj increases through 1, 
the synchronous periodic orbit loses its stability via saddle-node or pitchfork bifurcation. 
On the other hand, when Xj decreases through —1, it becomes unstable via period-doubling 
bifurcation. 

It follows from the coupling condition (|5.3|) that 


N 


= 0 . 


(5.17) 


1=1 


Hence the Lyapunov exponent cxo becomes 


1 


m—1 




m —^oo jji 


(5.18) 


t=0 


It is just the Lyapunov exponent of the uncoupled ID map. While there is no coupling effect 
on (Jo, coupling generally affects the other Lyapunov exponents (j/s (j ^ 0). 

In case of the global coupling of the form (|5.5|), the reduced coupling functions become: 


* \ V'(A') for I # 1, 


where 


V(X) = -r'(X). 


(5.19) 


(5.20) 


Substituting into Eq. ( b.l4| ), we obtain: 


m—1 


aj = 


lim — In \u'{Xl) \ for j = 0, 

T m—1 

- «''W)I tori 7^ 0- 


(5.21) 




t=o 


Here (Jq is just the tangential Lyapunov exponent cr|| of Eq. (|4.1CI|) , characterizing the mean 
exponential rate of divergence of nearby orbits along the synchronization line (|5.8|) . On 
the other hand, all the other ones aj {j ^ 0) (characterizing the mean exponential rate of 
divergence of nearby orbits across the synchronization line) are the same as the transversal 
Lyapunov exponent (j^ of Eq. (|4.11|) , i.e.. 


(Ji — ■ ■ ■ — ajv_i — cr_i_. 


(5.22) 
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Consequently, there exist only one independent transversal Lyapunov exponent a±. 

As for the case of two coupled maps in Sec. |I^, we choose u{X) = 1 — AX"^ as the 
uncoupled ID map, and consider a global-coupling case with r{X) = m(X), i.e., a global- 
coupling of the form 

c ^ 

t.(A'(l).....A'(iV)) = ^i:WX(m))-«(A'(l))] for]V>2. (5,23) 

[Here the case for = 2 is just the dissipative coupling of Eq. (^^) .] For this kind of global 
coupling, the transversal Lyapunov exponent ( b.22|) is given by 


a_L = cr|| -I- In 11 — c| 


(5.24) 


As in Sec. we also consider the intermittency associated with a saddle-node bifur¬ 
cation to a pair of synchronous periodic orbits with period 3. Then the phase diagram in 
Fig. g holds for all N {N > 2) globally-coupled maps, because the two independent Lya¬ 
punov exponents a\\ and are the same, irrespectively of N. Thus there exists a critical 
line segment joining two ends Q = 0 and = 2 on the synchronous saddle-node-bifurcation 
line A = Ac = 1.75 in the c — A plane for any N globally-coupled case. The critical be¬ 
haviors of N globally-coupled maps for N > 2 are also essentially the same as those of 
two dissipatively-coupled maps, in which case there are two kinds of critical behaviors (for 
details of the N = 2 case, refer to Sec. 


m 


C. Renormalization analysis of many globally-coupled maps 


As seen in Sec. [VB| , there exists the same critical line segment associated with inter¬ 
mittency for N [N = 2,3,4,...) globally-coupled maps, irrespectively of N. Following the 
same procedure of Sec. 0 for two coupled maps, we extend the renormalization results of 
two coupled maps to arbitrary N globally-coupled maps. We thus hnd two hxed maps of 
the reduced renormalization transformation and obtain their relevant eigenvalues. One hxed 
map with no relevant CE’s is associated with the critical behavior inside the critical line 
segment, while the other one with relevant CE’s is associated with the critical behavior at 
both ends. 

To study the intermittency in the vicinity of a saddle-node bifurcation to a pair of 
synchronous orbits with period p in an iV-coupled map M of Eq. ( 0 ), consider its pth 
iterate 


: Xt+i{m) = W^P\a^-^Xt) 

= W^P\X,{m),X,{m + 1),..., W(m - 1)), 

m = 1,. .., iV, 

where satishes a recurrence relation 


(5.25) 


W(P)(X) = W{W^P-^\X), W(^’-^)(ctX), ..., W(^’-^)(cT^-^X)), (5.26) 

and it can be also decomposed into two parts, the uncoupled part and the remaining 
coupling part, i.e.. 
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(5.27) 


H/(i.)(X) = !iW(x(l)) + |W'<r>(X) - !i<'>(A'(1)|. 

For the threshold value Ac of a synchronous saddle-node bifurcation, p synchronous hxed 
points of appear. Shifting the origin of coordinates (X(l),... ,X{N)) to one of the p 
hxed points (X*(l),.. .,X*{N)) (X*(l) = ■ ■ ■ = X*{N) = X* = u^p\X*) for A = Ac), we 
have 


where 


T : Xt+i(m) = F(cr'" ^x^)) 

= f{xt{m)) -h g{xt{m),..., Xt{m - 1)), 
m= 1,... ,N, 


f{x{l)) = u^P\x{l)+X*)-X*, 

g{x) = +X*,...,a;(iV) +X*) +X*). 


(5.28) 

(5.29) 

(5.30) 


Here x = {xi,... ,xn), the uncoupled part / for the critical case A = Ac satisfies the 
condition (|3.7| ), and the coupling function also obeys the condition (|5.3|). 

We employ the same renormalization transformation M of Eq. (|3.9|) with the rescaling 
operator al, where a is a rescaling factor, and I is the N x N identity matrix. Applying the 
renormalization operator Af to the X-coupled map ( b.28| ) n times, we obtain the n-times 
renormalizaed map T„ of the form. 


T„ : Xt+i{m) = F„(a”" ^x*) 

= fn{xt{m)) -h gn{xt{m ),..., Xt{m - 1)), 

m = 1,..., X. 


(6.31) 


Here the uncoupled and coupling parts / and g satisfy the following recurrence relations: 


/„+,(j;(l))=a/„(/„(hh)), (5,32) 

a 

iV-1 

(?.+i(x) = aUFni-)) + agr^iF^i-), ..., F^{- -^)) 

CK a CK 

-a/„(/„(—)). (5.33) 

a 


Then Eqs. (|5.32| ) and (|5.33|) define a renormalzation operator TZ of transforming a pair of 
functions (/,(?): 


f fn+l \ 

\ gn +1 j 


= n 


fn 

gn 


(5.34) 


A critical map Tc is attracted to a fixed map T* under iterations of the renormalization 
transformation M, 


T* : xt+i{m) = F*{a^-^^t) 

= f*{xt{m)) 

+g*{xt{m),xt{m + 1 ),..., Xt{m - 1)), 

m = l,...,X. (5.35) 
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Here {f*,g*) is a fixed point of the renormalization operator 71, i.e., = 7Z{f*,g*). 

Since f* is just the ID hxed map (|3.16|) , only the equation for the coupling hxed function 
g* is left to be solved. 

As in case of two coupled maps, we construct a tractable recurrence equation for a 
reduced coupling function of gix.) dehned by 


- Sr 


(5.36) 


x{l)=---=x{N)=x 


because it is not easy to directly solve the equation for the coupling hxed function. Differ¬ 
entiating the recurrence equation (|5.33|) for g with respect to x{l) {I = 1,..., N) and setting 
a;(l) = ■ ■ ■ = x{n) = x, we obtain: 




2 = 1 


a a 

N 


(5.37) 


N 


The reduced coupling functions G^^^’s satisfy the sum rule ( p.!?! ), i.e., '^G^’'\x) = 0, and 


1=1 


they also satisfy (x) = (x) due to the periodic boundary condition. 

As shown in Eq. (|5.19|) , there exists only one independent reduced coupling function 
G(x) for the globally-coupled case, such that 


G(2)(a;) = ... = G(^)(^)=G'(a^), 

= (l-iV)G(a:). 


(5.38) 


Then it is easy that the successive images {G((^(x)} of {G*^^^(a;)} under the transformation 
( f).37| ) also satisfy Eq. (|5.38| ) [i.e., G(f^ = •.. = G[fy^(x) = G„(x), G()^)(x) = (1 — A^)G„(x)]. 
Consequently, there remains only one recurrence equation for the independent reduced cou¬ 
pling function G(x) 


G.+i(fy = [/'(/.(-)) - iVG.(fy(-))]G.(-; 

OL a a 

+Gn{fn{-))a-). 

a a 


(5.39) 


Then, together with Eq. ( |5.32| ), Eq. (^.39|) dehnes a reduced renormalization operator 7Z of 
transforming a pair of functions (/, G): 


fn+l 

Gn+1 


= n 


fn 

Gn 


(5.40) 


Since the reduced renormalization transformation (|5.4CI| ) holds for any globally-coupled cases 
of > 2, it can be regarded as a generalized version of Eq. (|3.20|) for the case of two coupled 
maps. 
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A pair of critical functions (/c, G) is attracted to a pair of fixed functions (/*, G*) under 
iterations of 'JZ. Here f* is the ID fixed map ( p.l6| ), and G*{x) is the independent reduced 
coupling fixed-function of 5 '*(x) [i.e., G*^‘^'>{x) = ■■■ = G*^^\x) = G*{x), G*^^\x) = (1 — 
N)G*{x)]. As in the two-coupled case (iV = 2) of Eqs. (|3.21|) and (|3.22|) , we hnd two series 
solutions for G*: 


G*j{x) = ^[1 zax"" ^ + z{z - H-], 

1 37 - 

G*ii{.x) = —[hax^-^ + 


J! _ + .. .1 


(5.41) 

(5.42) 


Here a and h are arbitrary constants. The solutions for G* have a common factor and 
hence the function N G*{x) becomes the same, independently of iV; this can be also easily 
understood by looking at the structure of Eq. (|5.39|) . In case of G/(a;), we can sum the series 
and obtain a closed-form solution, 


G',(x) = 


(5.43) 


However, unfortunately we cannot sum the series in G*jj{x) except for the cases 6 = 0 and 
z where we obtain closed-form solutions. 


G*jj(x) = 


0 for 6 = 0, 
j^[f*\x)-l] ioTb = z. 


(5.44) 


Linearizing the renormalization transformation TZ at its hxed point {f*,G*), we obtain 
a reduced linearized operator C transforming a pair of inhnitesimal perturbations {h, <h): 


hn+l 

4 * 71+1 


= c 


hn 

4*7. 


Cl 

As 


0 

C 2 


4*7 


(5.45) 


where 


hn+i{x) = [Cihn\{x) = ar\f*{-))hn{-) + ahn{f*{-)), 

a a a 


^n+l(x) = [£24>n](x) [Ciihn]{x), 


.X 


. X 


X , 


.X , 


.X . 


.X, 


[C- 2 <in]{x)= [/•'(/•(-)) - AfG*(r(-))]4„(-) + [/•'(-) - lVG*(-)]4>„(r(-)), 

a a a a a a 

[c,K]{^)= [r\f\-))G\-)-NG*\f\-))G\-) + G*\r{-))r\-)]K{-) 

a a a a a a a 

+G*{-)K{n-))+G*{n-)m-). 

a a a a 


(5.46) 

(5.47) 

(5.48) 

(5.49) 


It follows from the reducibility of C into a semiblock form that to determine the eigenvalues 
of C it is sufficient to solve the eigenvalue problems for Ci and C 2 independently. Then the 
eigenvalues of both Ci and C 2 give the whole spectrum C. 

The eigenvalue equation for Ci is given by Eq. (p.31|) . As mentioned there, that is just 
the eigenvalue equation for the ID map, in which case the complete spectrum of eigenvalues 
and the corresponding eigenfunctions have been found in Refs. 0. 
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We next consider an infinitesimal coupling-perturbation of the form (0, <h) to a fixed 
point (/*,G*). If an independent reduced coupling perturbation d)* satisfies 


A$*(a;) = [£2<h*](a;) 


,x. 


.X, 


,x. 


= [r{r{-)) - NG*{r{-)w{-) 

a a a 

a a a 


(5.50) 


then it is called the independent reduced coupling eigenfunction with CE A. Note that the 
eigenvalue equation ( |5.5CI|) for any N becomes the same as that for the two-coupled case 
{N = 2), because the function NG*{x) is the same, irrespectively of N. Hence we follow 
the same procedure in Sec. 0 for two coupled maps, and hnd the same relevant CE’s for 
any N globally-coupled maps as follows (for more details, refer to Sec. |IITD : 


(1) G*(a:) = G}{x) 


There exist no relevant CE's. 

(2) G*{x) = G*j{x) 

There exist {z - 1) relevant CE's such that Xk = 

with eigenfunction ^l{x) of Eq. (3.37) {k = 0,z — 2). 


(5.51) 


(5.52) 


As in Sec. [VB| , we also consider the global coupling of Eq. ( b.23|) with r(X) = u{X) = 
1 — AX^ and study the intermittency associated with a saddle-node bifurcation to a pair 
of synchronous orbits with period p = 3. Considering the third iterate of M [see 
Eq. (|5.25|) ] and then shifting the origin of coordinates to one of the three synchronous hxed 
points for A = Ac, we obtain a map of the form T (|5.28|). The uncoupled part / has the 


form ([4.5|) , ans hence that corresponds to the most common case with the tangency-order 
z = 2. The independent reduced coupling function of the coupling part g{x.) is also given by 


G{x) = e = c^ -3c^ + 3c. 


(5.53) 


Consider a pair of initial functions (/c, G) on the synchronous saddle-node bifurcation 
line A = Ac, where fc(x) is just the ID critical map and G(x) = ■§ fc(x)- By successive 
applications of the reduced renormalization operator 7Z of Eq. (|5.4CI|) to (fc, G), we have 


fn(x) = afn-l(fn-l(-)), Gn(x) = ^ f^(x) 

a JXI 


(5.54) 

(5.55) 


where the rescaling factor for z = 2 is a = 2, /o(x) = fc{x), Gq{x) = G{x), and Cq = e. Here 
fn converge to the ID hxed function f*{x) of Eq. ( |3.16|) with z = 2 as n —>• cx). 


As shown in Sec. any initial e inside the open interval (0, 2) converges to the super¬ 
stable hxed point e* = 1 under successive iterations of the transformation ( |5.55|) . The left 
end of the interval is un unstable hxed point e* = 0, which is also the image of the right 
end e = 2 under the transformation ( |5.55|) ; all the other points outside the interval diverges 
to the minus infty under iterations of the transformation (|5.55|) . One can see easily that 
the interval [0, 2] of the parameter e corresopnds to a critical line segment joining two ends 
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Q = 0 and Cr = 2 on the synchronous saddle-node bifurcation line A = Ac in the c — A plane. 
Hence any initial G{x) inside the critical line segment is attracted to the first independent 
reduced coupling hxed function G*j{x) = jff*{x) under iterations of TZ, while G(a;)'s at both 
ends are attracted to the second independent reduced coupling fixed function G*jj{x) = 0 
with 6 = 0. 

VI. SUMMARY 

The critical behaviors for intermittency in two coupled ID maps are studied by a renor¬ 
malization method. We hud two hxed maps of the renormalization transformation. Although 
they have common relevant eigenvalues associated with scaling of the control parameter of 
the uncoupled ID map, their relevant CD’s associated with coupling perturbations vary 
depending on the hxed maps. These two hxed maps are also found to be associated with 
the critical behavior near a critical line segment. One hxed map with relevant CD’s govern 
the critical behavior at both ends of the critical line segment, while the other one with no 
relevant CD’s governs the critical behavior inside the critical line segment. We also extend 
the results of the two coupled ID maps to many globally-coupled ID maps. 
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FIGURES 


FIG. 1. Plots of the curve = 2en-i — and the straight line = Cn-i- 

FIG. 2. Phase diagram for a dissipative-coupling case (4.2). Here solid circles denote the data 
points on the (T_l = 0 curve. The region enclosed by the = 0 curve is divided into two parts 
dnoted by P and C. A synchronous period-3 (chaotic) attractor with cj|| < 0 ( cj || > 0) exists in the 
part P (C). The boundary curve denoted by a solid line between the P and C parts is just the 
critical line segment. 

FIG. 3. Iterates of the dissipatively-coupled map M of Eq. (2.2) with the uncoupled ID map 
(4.1) and the coupling (4.2) for (= Ac — A) = 10“^ and c = 0.5; (a) an orbit with initial point 
(0.1, 0.3) is attracted to a synchronous chaotic attractor on the synchronization line Y = X after 
exhibiting a short transient behavior, and (b) the plot of the A-component A) ' of the third iterate 
of M versus a discrete time t shows intermittent alternations between laminar and turbulent 
behaviors. 

FIG. 4. Plot of (t_l (= In |1 — c|) versus c for A = Ac- The values of a± at both ends of the 
critical line segment are zero, which are denoted by solid circles. 

FIG. 5. Transition from a synchronous to an asynchronous state for A = Ac. A synchronous 
period-3 orbit denoted by circles is no longer an attractor outside the critical line segment, and 
new asynchronous attractors exist for (a) c = —0.0001 and (b) c = 2.0001. We denote a pair 
of asynchronous period-3 attractors by uptriangles and downtriangles in case of (a), while the 
asynchronous attrractor with period 6 is denoted by “stars” in case of (b). 
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